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Morphological, karyological and molecular delimitation 
of two gentians: Gentiana crassicaulis versus G tibetica 
(Gentianaceae) 
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Abstract Gentiana crassicaulis Duthie ex Burk. and G. fibetica King ex Hook. f. are 
confusing gentians. The former has smaller corolla and blue-purple limbs whereas the latter 
has larger corolla and purple-brown limbs. Intermediate types are often found in southeast 
Xizang where both species are putatively considered to occur. A question has been raised 
whether these intermediate types were produced by introgression between the two species, 
which further contribute to the taxonomic ambiguity. In the present study, we comprehensively 
compared morphological variations, chromosome numbers, and sequences of two DNA 
fragments from populations of the intermediate types and those of the two typical species. We 
sequenced cpDNA tmS-G and nuclear ITS fragments of 69 individuals of 10 populations, 
counted chromosome number of 26 individuals of six populations and measured 
morphological variation of 54 specimens. We found that corolla size is highly correlated with 
ploidy level. The populations of typical G crassicaulis from Yunnan were diploid (2n=26) 
with small corolla while those of G tibetica from western Xizang were tetraploid (2n=52) 
with large corolla. The populations of the intermediate types from southeast Xizang were 
tetraploid with large corolla, and a few individuals with purple to brown limbs. These features 
suggest that the intermediate types show more similarities to G. tibetica. However, DNA 
sequence evidence unambiguously grouped those populations with G crassicaulis. Both 
nuclear and chloroplast DNA showed consistent results: the intermediate populations showed 
identical cpDNA and very similar ITS sequences to typical G crassicaulis. No sequence 
variation beyond what occurred in typical G crassicaulis and G tibetica was found from 
multiple individuals of the intermediate types. Thus the hypothesis of a recent introgression or 
hybridization between the two species could be rejected. Our results demonstrated that 
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polyploidization had contributed to the intraspecific variation in G crassicaulis and that the 
tetraploid individuals of this species showed morphological convergence with tetraploid G 
tibetica, which blurred the morphological delimitations of the two species in southeast 
Xizang. Nevertheless, molecular markers can discriminate convergence and correctly attribute 
the intermediate populations or individuals to the correct species circumscription. Our results 
provided a case study to circumscribe the closely related species based on a comprehensive 
combination of morphological, karyological and molecular evidence. In addition, the 
geographical distributions of diploids and tetraploids within G crassicaulis and the 
implications for the origin and dispersals of this alpine species were also discussed. 

Key words Gentiana crassicaulis, Gentiana tibetica, polyploidization, morphological 
convergence, species delimitation. 
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O1 םם‎ [] 0 0 0 D 
Table 1 The origin of the materials, voucher information and the number of individuals used for DNA sequencing 
D D D D 9 Ek? 
Population Locality Latitude Longitude Altitude Voucher” Sample 
code (m) size” 
LT [] 30°09' E 99°14'N 3800 ה‎ , [| [| D] CY. M. Yuan 5 
Litang, Sichuan & X. J. Ge) CN2K2-92 (IBSC) 
LJ [] 26°53'E 10015'እ 3200 |] (J. Q. Liu) 2286 5 
Lijiang, Yunnan (HNWP) 
RM [] 29°39' E 98°27' N 3700 [| , [] (Y. M. Yuan 10 
Rumei, Markam, & X. J. Ge) CN2K2-126 
Xizang (IBSC) 
BD [] 208 98°41'N 3400 ole [] (Y. M. Yuan 10 
Bamda, Markam, & X. J. Ge) 2002-001 (IBSC) 
Xizang 
MK [] 29°42'E 98°36'N 3800 ₪ [] (Y. M. Yuan 10 
Markam, Xizang & X. J. Ge) 2002-002 (IBSC) 
ZG [] 29°55'E 97°39' N 3820 [| , [] (Y. M. Yuan 10 
Zogang, Xizang & X. J. Ge) CN2K2-122 
(IBSC) 
BM [] 29°57' E 94°54’ እ 3075 |] J. Q. Liu) 1214 B 4 
Bomi, Xizang (HNWP) 
NG [] 29°03' E 90°24" N 3700 [] (.. Q. Liu) 1117 5 
Nagarzé, Xizang (HNWP) 
Qx [] 29°18' E 90°21'N 3650 |] (J. Q. Liu) 1110 5 
Qüxü, Xizang (HNWP) 
QS [] 29°02' E 92*18' እ 3980 |] (J. Q. Liu) 1156 5 
Qusum, Xizang (HNWP) 
11 Du [] D aBSOJ D D [] [] በ] (HNWP); 2) 
0 [] םםםםםםםםםה‎ D OO 
1) Voucher specimens were kept in South China Botanical Garden herbarium (IBSC) and Northwest Plateau Institute of 
Biology (HNWP), the Chinese Academy of Sciences. 2) The numbers denote the individuals being sequenced in each 
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በ] [] በ 72°C7minf በ በ | በ በ D በ በ CASpure PCR[T D] CU C] D (1 D] D. U | D OH 
םםםםהםםםהםםהםםםםהםםיםםםםםםםםם‎ በ በ Amersham] UO 
DYEnamic Dye Terminator Cycle Sequencing Kit[] [] [] ה‎ 1 [| [| E] ቨ] [] [] D]. השתהה‎ 
מםםםם.ם‎ በ በ Megabase 5000 | [| U D [በ በ በ በ DNAD | D. ITSE] D] D በ U በ 
GenBank] [O O [] [] [] G frigida Haenke[] ITSO O [] (Yuan et al., 1996)[] O O trnS-G[] O 
OO [] [] Glycine dolichocarpa Tateishi & Ohashi] O O O [| [| [| (Sakai et al., 2003)[] [O 
םם:.םםםםםםטםםםםםםםםםםםםם‎ 0 0 GenBank ([] 3)[] 
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The geographical distribution of the populations studied in Gentiana crassicaulis and G. tibetica and their 








O00 Gerassicaulis, 2n=2x/4x=26/52 
OOD G.tibetica, 2n=4x=52 
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በ1 


Fig. 1. 


corresponding chromosomal ploidy. Population codes referred to Table 1. 
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[ [ [ ו‎ 0 
Table 2 The populations studied, individual numbers and their chromosome numbers 
oo” 0 ו‎ 
Population” The number of analyzed individuals Chromosome number 
LJ 1 2n = 26 (Yuan et al., 1998) 
RM 10 2n= 26 
BD 9 2n = 6 
እብሩ 1 28 = 52 
ZG 2 2n=52 
BM 2 2n=52 
LK 2 2n = 52 
L 10 












































































































































































































































DO 


1) Population codes referred to Table 1. 








03 00 [] ITS] trns-G| D [] 
Table3 The information of ITS and trnS-G sequences of G. crassicaulis and G. tibetica 
[] DNA / 
fragment Species GenBank Length (bp) Substitution Indels (bp) Genetic distance 
accessions P? (96) 
: G. crassicaulis DQ398636 626 
nrDNA: ITS crassicaulis Q 9 10) 10 
G. tibetica DQ398640 629 
cpDNA: trnS-G G. crassicaulis DQ398743 528 
EE 3 4 (36) 0.6 
G. tibetica DQ398747 567 
POOP በ ה‎ 000080 


























11061 was regarded as a nucleotide substitution. 




































































CLUSTAL XQ O (Thompson et al., 1997) | | || | | ||] ] 1 []‏ |[ ה ם ה ם 
D D UO D U . 0 1 ሄጩ)[ [1 0 0 D] በለ ቓቓሠጵ/ጩመበበበበበበበበገበበበ[0በ00000በ‏ 
በ PAUP4.0b10(Swofford, 2000)] 0 ] D] D] C] በ D] በ በ በ (MP)‏ .םםםםםםםםם 
ቨ ቨ ፡‏ 1 0 1 1 1 ][0 0 0 1 0 1 1 1 1 1 ה ][ ][ |[ |[ |[ [ םה OOOO ANDO‏ 


P* (54): Uncorrected p-distance, an 
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MPU በ] [| ፲ | ፲ D  ፲8፻ቨ D] ገ በ ል MULTARS, ACCTRAND ቨ , 1000] | HE | ገ በ U UI 
םםםםםםםםםם‎ ሀ በ (bootstrap, םאה‎ በ 11111] በ [ በ [ 
0“ P” OOO 00 በ MEGA3.0(Kumar et al., 20040 000000000000000 
በበበ በ פםםםםםםםםםםםםם‎ 


2 DO 


21 00000 

በ1161: 1 1 20:00በ:00በ0በ0.000ሀ0000›፡ מםםםםםטמםם‎ 
(2n=2x=26), 70 OOO [| 52, በ] 1 1 D 21 un BUD HD DH BH D HD DO 0. 
םםטםםםםםםהםהלהםםםםםםםםםם‎ በ ሀ ሀ ሀ ሀ (GMJL] [ [| BD) 
በ ሀገ በ Mk), 00 20; [ ቨ , 0 በ ፳እብ] .םםםםםם הרת‎ በ በ፲ በ በ ሀ ሀ ቨ 
םםםםםםם‎ ዐቭበበቨበበበበበቨበቨ በበበ ሀን በበ በ፡በ[;ቧበ ;ቨ [ከ ቨ 
በ 0 በ ה האוא‎ BMO 50 [| DD ገ [| ፲ በ Crs. 000000000000000., 
םםםםםםםםם‎ 

22 סםםםםםה‎ 
הלמםםם.םםםםםהםםםהםהםםהםםםםםםהאהםםםםםםםםם‎ 
םםםםםםםםםםםםנםם‎ 2.0-2.5 c0 DD. ሀ ሀ በ «Mit ceu uu 
הפא‎ 51 1111110111 ፲ በ ሀ 
በ], 38-45 em] D. 110100: 
በሀ በ [፲ ቨ ፲ ሀ በ D בש‎ 
በ1000: በ10[በ10[0ህቭ. 
8.0-19.5 cml] „00000000 
סםםם םםםם‎ 
פםםםםםםםםם‎ 
0300.00 YY [| (SOU [] 
באו םמםםםם‎ D B D OU UI 
בבםםםםהםהםהםהפאםםםם‎ 
በሀበ በ፲ ሀ ם‎ EMU ב ם‎ 
(:)1 1 ששה‎ በ ][1 | ገ | ; 
በ GOLD [| ፲ U (እብን[ DU U ፲ ፲ U [፲; 
הםםםםםםםםםםם‎ 00900 
פםםסםםם םםםםםםםםםם‎ 
פםםםםםםםםםםםםטםםםם‎ 
በበ] ፲ገበበበ በሀ; በበ 
םםOםOOOOOOOOOOOOGO‎ 
סםםםםםםםםםםטםם .םמ‎ 
በበበ פםםםםםםםםםםםם.‎ 
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Fig. 2. The scatter plots of corolla length against the 
inflorescence bract of individuals with chromosome 
counted. Solid triangles, typical G. crassicaulis and 
those collected from southeast Tibet with 2n=26; solid 
rhombi, typical G. tibetica; and blank triangles, those 
intermediate individuals collected from southeast Tibet 
with 2n=52. 
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20 מםםםםהםםםהםםםםהםםםםם‎ 
v " ממםםםם:םםםםםםםםםםם‎ 
TD מםםםםהםםםםהםםםםםםםם‎ 
«i | ine በ00000000000000003 
EE a ልቀ מםםםםהםםםהםהםםםםםםםם‎ 
me? ol Be, ,מםםםםםהםםםםהםםםםם‎ 
r3 5 מםםםםהםםםםםםםהםםםםם‎ 
4 000., [[ [ [ [ 0 2.6-3.5 emf OU 
| ממםםם:;:םםהםםםהםםםםםם‎ 
4 מממםםטם;תםםההמהת = =< א‎ 
CE ,םםםםםםםםםםםםםםםםם‎ 
מםםםםםםםםםהםםםםםםם‎ 
pac ን ንእ ን - , םמםםםםםםםםםהםםםםםםהם‎ 
D o0000000C ng 23 DNADDDU 
oe Tr CAA Ge םםםםםםםםםםםםםםם‎ 
109000 በ[ T nn 1 በ] 1100 1ቨ eed 0 በ በ [፲ሎቨ [ o 
a 3. The iones nus ERN length cee = in : 7 0 0 : : : : SE 
inflorescence ract o identifiec peice = a - , 4p" 
seg ל‎ iberica m southeast Tibet 00300400 00000000010 
wih Sess color records and אש אא‎ © 10. 
crassicaulis without fresh color records. [] [] [] [] [] [] ה‎ LH [] LH [] [] ITSO [] [] LH 
[] O 626 bp[] 629 bp, [] |] በ D] [ 100, 
OO TSiO 1፲521 םהה‎ 01 101,0በ0001[00000[ [0 [በ በ 1.05:1[1] 0000000 
በበ] በበበ በበገበበመጩጠበበበበበበበበመጠበበበበበ[በ[;በበበበ[[ D 
በበ] በ በበበበበ በ በበበበ በበ ገበበበበበ[[በ በበ [ በ msi [ ሀ[] [በ [ [ 
በ4 Dn ITS] ¢rnS-GQ [] D OU 
Table 4 The divergent sites of ITS and trnS-G sequences between G. crassicaulis and G. tibetica 
ITS 
1 1 4 5 5 5 6 
2 4 1 8 2 5 6 8 0 
9 2 5 7 4 3 7 2 5 5 
G. crassicaulis G 0 0 ፐ 0 A A T T 
G. tibetica A A T + G T G G C C 
2 2 4 כ‎ 5 
6 2 9 5 8 3 3 
8 9 3 8 6 3 4 
G. crassicaulis ה‎ ፐ 0 T A 
G. tibetica ል T G G [] ር: 6 


ቆ|] TGA, A] ATAGATATTATATGGATATGTACAC,[] ך‎ GAGGATACT. 
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מםםםםםםהםםםםהםםםםםהםםהםםםםהם 
םםםםםםםהםםהםהםהםם הש በበበ በ በበበ በ በ09 00 በ በ በ‏ 


Uu‏ פםםםםםםםםםםם 


LU 





[| חנ‎ [| [| [] (D 4), CE1.0(CF1.0, excluding uninformative 


CLUSTAL XJ [] [] በ D] 563 ሠር] D D D D] 
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ሀ 15;‏ םםםםםםםםםםםםטםםם ההםםהםהםאתהפפטםםםם 
הםםםם םםטםםםםםםםםםםטםםםםטםםטםםםםםםםם 
300000 00 | הפו bpi 567 bp, O0 |] [| LU] LI DI 0 ብ፲ [| 36 bpi OU‏ 528 


(substitution)([] 4), [] 
בםםםםםםםםם‎ 





מםםםםהה 24 


חב ITSO‏ || ה םה הא CLUSTAL‏ ה 





10191, 0 [1 ሀ በ ቨ 


characters), RI-1.0, [] O =18[[ tmS-G] O O 


በ ሀ በ በ polyAQ polyTO እ በገ በ በ በ በበበ ,ሀ110 1121, [0 ቨ በ ሀ ሀ DU 
10 0 [| ፲ @ 5), ርሃ=1.0, R7-1.0, [] ፲ =20[ ם‎ ITSD םםםם‎ D ፲ «ms-GD HE. D U D U 0[ 
םםםםםםםםםםםםםם‎ 00020. 0000000 ሀ ሀ ሀ GATS. 970 98; 
trnS-G, 88[| 98), 1] 1111111101000100ሀ0ዐ D U DD] U] D. HL D. ሀ 0 0 LILY 
הפפחאא‎ 10101110111,10,:0:80በ00[0ሀ000090000003ቨ ሀ በ ሀ ሀ ቨ 
האוא‎ ZG] BM] ITSO DO ወ ጠበ UH BOO DO ሀ በዐ ሀ ሀ በ ሀ ሀ በ በ OSD ህጻ 
NG 30000000000 RMQ BDOOOOOTSSOOOOOOOOOOOOOO 
በበበ מםםםםמםםםםםשםםםםםםםםםםטםטקםםםםםםםםם.‎ 
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Fig. 4. The single maximum parsimonious (MP) tree based on ITS sequences. Numbers above internodes are bootstrap 
values based on 1000 replicates. 
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Fig. 5. The single MP tree based on trnS-G sequences. Numbers above internodes are bootstrap values based on 1000 


he numbers in the brackets are the chromosome numbers of the corresponding populations analyzed. 





replicates. T 


300 
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OOOOGOCGOOOOOOOOOOOOO.OOOOOOOOG 
LU] D U B D OO (Grant, 1963; Arnold, 1997 00.0000000 
UOoO.OGOOOOOOOOOOOOOOOOOOOOOOOGO 
םםםםםםםםםםםםםםםםם‎ ቨ በ [ሀ LN L (Amold, 1997) 
םם.םםםםהםםםהםםםהםםםםםםםם םםםם‎ 
םםםםםםםהםם הםםטטםםםםםםם‎ mil. 


L3 
= 


- 
LJ E3 E AO NAA ፻፲ Ilt. 





OU 
םםבצשםםםםםטםםםםםםםהםםםםםםםםםםם‎ 
ם‎ (Mogensen, 1996), 1 1 000000 
םםםהפעם:;םםםםםםםםהםהםםםםםםםםםםםם‎ 
םםםםםםםםםםםםםםם םםםםםםםםםםםם‎ 
[|] EL E] E] [| (additivity) (Wendel et al., 1995; Koch et al., 2003) 
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በገ በገበበ በበበ በበበሀንበሀበ በን በበ በበ በበበ በ[9በ[,;ቧ› ሀ 
በገ  በገበበበቨበ በበ በሀ;።ን።ሀበበበቨበሀበበሀ ገ በበ በበበ 
በገበበገበ በ በ በ በ 10ሀ0በ00በበ80ፌ0[800ሀ81000800ሀ0000090ዐቨ 
በገ በገ በ በበ በበበ፡ዐሀን በሀበ በበር ቧበበ[፲9በ0ሀቨ። ሀ ሀ ቨ 
በገ] በበበ በ በበ፡ሀበ በበበ ה‎ ODDUN 
פםםםםםםםםםםםםםםםםםם‎ 
םםםםםהםםםםםםםם‎ 00i 
פםםםםםםםםםםםםםםםםם‎ 
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םםםםםםםהםםםהםםםטםםםם 
פםםםםםםםםםםםםםםםםם 
U DU BU D B [በ U በ (Tautz et al., 2003), []‏ םםםםםםםםםםםםםםתםם 
በበ 1101001000 ሀ ሀ ሀ ሀ ሀ DD.  (erypti species) (Dick et al, 200‏ 



















































































- 
- 
- 
- 
- 
Ez 
L1 
E 
L1 
L1 
E 
L1 
L1 
L1 
E 
L1 
L1 
መ ET, ESTY ES ₪: EST FFT EI 
፲=1. ET E-T- EZT =. ב‎ ET: E-T- ESL TE)" E21 





0 

0 

[] 

[] ; 
በገ በ [በ በ םםםםםםםםם‎ 
[] 

[] 

[] 

[] 





(ታን 


































































































Dawson et al., 20050 (| [| [| [| 1 1 [| | D.U | | d || 1 | (11 1 D U D [ ን[| [| 
በ10[0[0[0[00[000000000000000000000000000[0[00[0 
በገበበበበገበበበበበዘበበበበሀበበበበበበበበበበበበበበ;በበ[በ 
(Yang et al., 2003; Luo et al., 2005; Wang et al., 2005; Liu et al., 200000 [| [| [| E] חח‎ E] D] 
םםםם‎ D] D] D] (Yuan. & Kipfer, 1997); םםםםםםםםםההםההםהםהםהםהםהםםה‎ 
በበ] .11[10 0451711011 10.110000በ010ሀ1000000000000ቨ 
[| [] [] [| [| [| (Zhang etal., 2006, unpublished data)[][] [| |] | | || | | |||] 111] [| [| 
በ በበ0በ0በ0በ0በ00በ000000000000000.0[01[0[01001[010[ቨ[ሸ[ [ ሀ] በ 
(EL OL EL OL ET CE DE CE D] UI D]. הפע‎ በ O cpDNA ם ב6-?‎ UU U ם ם ם‎ 154] 0.679] E] 
በ በ0በበበበበበበበበበበበበበበበበበበበበ።በበበበበበበበበበበበ[በ 
በበበበ[00በ0000000000000000000000000000000ቨ 
በበበ በ በበበበበበበበበበበበበበበቨበበበበበ።በበበበበ[በበበበ[በ 
በበበበበበ በ 11]101፲[000በ0በ0በበበ0በበበበበሀበበበበበበበ0በ0000በ00› 08 
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